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I . INTRODUCTION 


The  Hugoniot  of  rolled  homogeneous  armor  (RHA)  has  been  measured  as 
part  of  a coordinated  effort  to  determine  the  dynamic  and  quasi-static 
properties  of  steels  which  are  used  for  armors  and  penetrators.  This 
study  of  RHA  has  benefited  from  preceding  shock-wave  studies  of  iron. 
Barker's  reinterpretation  of  the  alpha-phase  Hugoniot  of  iron1  and  his 
collaboration  with  Hollenbach  in  the  shock-wave  investigation  of  the 
alpha-epsilon  phase  transition  in  iron2  have  been  invaluable,  providing 
experimental  techniques,  analytical  procedures,  data  for  comparisons,  and 
interpretations.  Further  information  about  the  shock- induced  alpha- 
epsilon  phase  transition  in  iron  has  recently  come  from  a careful  experi- 
mental investigation  by  Forbes.3 

RHA,  like  iron,  transforms  from  the  low-pressure  alpha  (bcc)  phase 
to  the  high-pressure  epsilon  (hep)  phase  at  a shock  stress  close  to  13 
GPa.  Just  below  the  transition  stress  impact  loading  produces  a two-wave 
stress  profile  which  consists  of  an  elastic  wave  followed  by  a slower 
Plastic  1 (PI)  wave.  Just  above  the  transition  stress  the  PI  wave  is 
followed  by  a Plastic  2 (P2)  wave  which  is  associated  with  the  epsilon 
phase.  The  stress  profile  reverts  to  a two-wave  structure  at  higher 
shock  stresses  where  the  velocity  of  the  P2  wave  exceeds  the  maximum  PI 
value.  Finally,  a single  shock  wave  is  propagated  at  higher  stresses 
where  the  velocity  of  the  P2  wave  exceeds  the  velocity  of  the  elastic 
wave. 


Hugoniot  data  for  the  alpha  phase  of  RHA  were  obtained  earlier  and 
have  been  reported.4  Additional  experiments  have  been  performed  to 
establish  the  Hugoniot  for  the  epsilon  phase  of  RHA.  This  report 
describes  the  experiments  and  analysis,  and  presents  data  for  the 
epsilon  phase  from  13  to  135  GPa. 


lL.  M.  Barker,  "Alpha-Phase  Hugoniot  of  Iron",  J.  Appl.  Phus. , Vol.  46, 
No.  6,  June  1976 , pp.  2544-7. 

2L.  M.  Barker  and  R.  E.  Hollenbach,  "Shock-Wave  study  of  the  ate 
Phase  Transition  in  Iron",  J,  Appl.  Phys.,  Vol.  45,  No.  11,  November 
1974,  pp.  4872-87. 

V.  W.  Forbes,  "Experimental  Investigation  of  the  Kinetics  of  the  Shock- 
Induced  Alpha  to  Epsilon  Phase  Transformation  in  Armco  Iron",  Naval 
Surface  Weapons  Center/White  Oak  Laboratory  Technical  Report  No.  77-137, 
December  1977. 

46.  E.  Hauver , "The  Alpha-Phase  Hugoniot  of  Rolled  Homogeneous  Armor", 
Ballistic  Research  Laboratory  Memorandum  Report  No.  2651,  August  1976, 
AD  #B012871L. 
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II.  MEASURED  PARAMETERS  AND  EXPERIMENTATION 
A.  Hugoniot  Equations 

The  thermodynamic  states  that  can  be  reached  by  shock  compression 
define  a curve  in  the  stress-volume  plane  called  the  Hugoniot.  In  the 
case  of  a'  multiple-wave  stress  profile,  the  stress  on  the  Hugoniot  is 
the  sum  of  the  stresses  carried  by  the  different  waves.  The  stress 
increase,  o-o^,  carried  by  each  wave  is  given  by  the  relationship, 

a-oi  - p^Cu-up,  (1) 

where  is  the  density  of  material  into  which  the  wave  advances,  U is 

the  shock  velocity  with  respect  to  material  ahead  of  the  wave,  and 
(u-ui)  is  the  change  in  particle  velocity.  The  specific  volume  on  the 

Hugoniot,  V,  is  given  by  the  relationship, 

V = V.  [l  - (u-u.)/u]  , (2) 

where  V^l/p^  In  the  case  of  iron  or  RHA,  the  stress  may  be  carried  by 
three  waves.  The  stress  in  the  elastic  wave  is 


poUeV 


the  maximum  stress  in  the  first  plastic  wave  is 

o,=a  + p U , (u  , - u ) , 
pi  e e pi  1 pi  eJ ’ 


(3) 

(4) 


and  the  maximum  stress  in  the  second  plastic  wave  is 


°p2  * V * v Up2  ("P2  - Upl>  • 


where  subscripts  o,  e,  pi,  and  p2  refer  to  unstressed  material,  the 
elastic  wave,  the  Plastic  1 wave,  and  the  Plastic  2 wave,  respectively. 
If  the  wave  velocity  is  measured  in  laboratory  coordinates,  the  expres- 
sion for  stress  in  the  ”P2  wave  becomes. 


ap2  ” °pl  + Ppl  ^p2 


Upl)tup2  - upl> 


The  epsilon-phase  Hugoniot  of  RHA  was  established  by  measuring  impact, 
shock  wave,  and  free-surface  velocities. 

B.  Test  Specimens  of  RHA 


Test  specimens  prepared  for  measurements  to  determine  the  epsilon- 
phase  Hugoniot  of  RHA  were  similar  to  test  specimens  previously  prepared 
for  measurements  to  determine  the  alpha-phase  Hugoniot.  Chemical 
analyses,  hardnesses,  and  preparation  procedures  reported  in  Reference  4 

8 


also  apply  to  specimens  used  in  the  current  tests.  Specimens  for  all  but 
one  of  the  current  tests  were  prepared  from  13-mm  thick  plate  stock;  for 
the  test  at  16.0  GPa,  the  specimen  was  prepared  from  100-mm  thick  plate 


3 3 

stock.  The  density  of  RHA  plate  stock  was  7.84  x 10  kg/m  . 

C . Experimental  Techniques 

Both  interferometric  and  photographic  techniques  were  used  for 
measurements  which  established  the  epsilon-phase  Hugoniot  of  RHA.  Peak 
stresses  up  to  34.5  GPa  were  obtained  by  plane  impact  experiments 
performed  in  a light-gas  gun.  Symmetrical,  RHA  -*■  RHA,  impacts  were  used 
at  peak  stresses  up  to  23.6  GPa,  but  a tungsten  impactor  was  used  to 
achieve  a peak  stress  of  34.5  GPa.  In  these  experiments,  the  free- 
surface  velocity  of  impacted  RHA  specimens  was  measured  continuously  by 
VISAR-type  instrumentation  (Velocity  Interferometer  System  for  Any 
Reflector)5,  and  a sequence  of  electrical  contactors  provided  a measure- 
ment of  the  impactor  velocity.  Explosively  accelerated  impactors  were 
used  to  obtain  higher  peak  stresses  of  63.6  and  135  GPa.  The  two-wave 
region  was  bypassed  in  going  from  a three-wave  shock  structure  at  34.5 
GPa  to  a single  shock  wave  at  63.6  GPa.  In  the  experiments  with  explo- 
sively accelerated  plates,  streak  camera  measurements  utilized  changes 
in  reflectance  to  detect  the  arrival  of  a strong  shock  wave  at  a free 
surface6.  Streak  camera  measurements  provided  the  shock  wave  velocity 
through  RHA  specimens  and  the  average  free-surface  velocity  of  a Type  304 
stainless  steel  plate  which  held  them. 


III.  ANALYSIS  AND  RESULTS 

Figure  1 shows  the  free-surface  velocity  profile  measured  in  the 
symmetrical  impact  experiment  at  a peak  stress  of  23.6  GPa.  The  elastic, 
PI,  and  P2  waves  are  identified  along  with  perturbations  believed  to 
result  from  elastic  reverberations  between  the  free  surface  and  approach- 
ing plastic  wave  fronts.  The  midpoint  of  each  wave  front  was  used  to 
establish  wave  arrival  at  the  free  surface.  Midpoints  are  indicated  by 
small  circles, and  times  t^,  t 2,  and  t4  identify  the  arrivals  of  the 

elastic  wave,  PI  wave,  and  residue  of  the  P2  wave,  respectively. 


Analysis  applicable  to  the  alpha-phase  Hugoniot  measurements  is 
included  in  Reference  4.  Briefly,  the  free-surface  velocity  profile, 
specimen  dimensions,  and  prior  knowledge  of  the  elastic-wave  velocity 


are  used  to  determine  the  velocity,  U , of  the  PI  wave.  For  the  elastic 


5L.  M.  Barker  and  R.  E.  Hollenbaoh , "Laser  Interferometer  for  Measuring 
High  Velocities  of  any  Reflecting  Surface ",  J.  Appl.  Fhys. , Vol.  43, 
No.  11,  November  1972,  pp.  4669-75. 

6G.  E.  Hauver  and  A.  Melani , "The  Hugoniot  of  5083  Aluminum ",  Ballistic 
Research  Laboratory  Memorandum  Report  No.  2345,  December  1973. 
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Figure  1.  Free-surrace  velocity  profile  measured  in  the 
impact  experiment  at  23.6  GPa. 
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wave  and  the  PI  wave  at  the  phase  transition,  the  change  in  free-surface 
velocity  is  assumed  to  be  twice  the  change  in  particle  velocity,  i.e.. 


Cuf  - ufi)  = 2(u  - Uj),  (7) 

where  subscript,  f,  refers  to  the  free  surface  and  subscript,  i,  refers 
to  the  initial  condition  (before  wave  arrival) . The  adequacy  of  this 
assumption  was  supported  by  the  investigation  of  alpha-phase  RHA  reported 
in  Reference  4.  Forbes3  cites  the  work  of  Walsh,  et.al.5  , and  Barker 
and  Hollenbach2  in  defending  the  adequacy  of  the  "factor-of-two"  assump- 
tion for  the  case  of  iron. 

The  epsilon-phase  Hugoniot  of  RHA  was  established  by  evaluating 
up2  and  Up2-  In  the  symmetrical  impact  experiments,  up2  = uq/2,  where  uq 

is  the  velocity  of  the  RHA  impactor.  In  the  experiment  with  a tungsten 
impactor  and  in  the  experiments  with  explosively  accelerated  plates,  up2 

was  determined  by  an  impedance  solution5.  The  shock  velocity,  U was 

P* 

determined  from  measurements  of  free-surface  velocity.  Assumptions 
involved  in  the  determination  of  Up2  are  explained  by  reference  to  Figure 

2,  a distance-time  plot  in  laboratory  coordinates.  Point  (X^,  t^)  must 

be  located  to  establish  Up2  (See  Equation  6) , where 


Up2  = ^X3  " X0^^t3  ' t(P 


It  can  be  shown  that. 


U2<X2  - X0>  * VX4  - X0>  - U1  U2ft4  * *2> 

V " u2(t4  - t())  . Uj(t2  - t())  - (x4  - k2) 

Time,  tQ,  is  established  from  the  specimen  thickness,  (Xj  - X^) , and 

prior  knowledge  of  the  elastic-wave  velocity,  U . In  the  expression  for 
* e 


r 

(X2  - X0)  - (Xj  - Xq)  ♦ J uf  dt, 


(x4  - X2)  = J uf  dt. 


(X4  - Xq)  = (X2  - Xq)  ♦ (X4  - x2). 


4 O'  v 2 
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Unlike  Forbes1 * 3 *  who  assumed  that  U^  and  U2  depend  on  the  elastic-wave 

velocity,  it  is  assumed  here  that  elastic  reverberations  can  usually  be 
ignored,  with 


When  the  residue  of  the  P2  wave  is  less  than  13  GPa 


where 


(From  Reference  4,  CQ  = 4510  m/s  and  S = 1.430).  Figure  1 offers  some 

basis  for  ignoring  the  elastic  reverberations  which  appear  in  each 
instance  as  a perturbation  below  the  midpoint  of  the  wavefront  where  the 
arrival  time  is  evaluated.  With  respect  to  material  into  which  the  P2 
wave  advances. 


Figure  3 shows  the  free-surface  velocity  profile  obtained  by  veloc- 
ity interferometry  from  the  experiment  at  34.5  GPa  in  which  a tungsten 
impactor  was  employed.  The  shock  velocity,  Up2,  was  again  determined  by 

Equation  8.  The  particle  velocity,  Up2,  was  located  at  the  intersection 

of  the  (ppj  Up2)  line  for  RHA  with  the  known  Hugoniot  for  tungsten7  which 

passes  through  zero  stress  and  the  impact  velocity,  uQ,  as  shown  in 

Figure  4.  The  (p^  U^2)  line  comes  from  Equation  5 and  in  this  case  is 

the  line  with  slope  (p  . U _)  which  passes  through  the  transition  state 


Epsilon-phase  Hugoniot  points  were  also  obtained  by  impedance  solu- 
tions in  the  tests  with  explosively  accelerated  plates.  The  measured 
shock  velocity  through  RHA  specimens  on  the  304  stainless  steel  plate  was 
used  to  establish  the  line  with  slope  (pQ  U^2)  in  the  o-u  plane  (See 

Figure  5) . This  line  intersected  the  release  adiabat  of  304  stainless 
steel  at  the  Hugoniot  state  in  RHA.  The  release  adiabat  of  the  stain- 
less steel  was  located  by  a measurement  of  free-surface  velocity,  and 
was  generated  by  assuming  a constant  yp,  where  y is  the  Gruneisen 


1R.  G.  McQueen,  S.  P.  Marsh,  W.  J.-  Carter,  J.  N.  Fritz , and  J, 

"The  Equation  of  State  of  Solids  from  Shock  Wave  Studies",  R 

ity  Impact  Phenomena,  edited  by  R.  Kinslow,  Academic  Press,  i 

1970,  Chapter  VII. 
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Figure  3.  Free-surface  velocity  profile  measured  in  the 
impact  experiment  at  34.5  GPa. 
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Figure  4.  Impedance  solution  for  the  impact  experiment  at  34.5  GPa. 
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Figure  5.  Impedance  solution  for  the  impact  experiments  with 
explosively  accelerated  impactors. 
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coefficient  and  p is  the  density.  The  data  used  for  304  stainless  steel7 
were:  pQ  = 7.896  x 103  kg/m3;  yQ  - 2.17;  U = 4569  + 1.490  u. 

Data  for  the  epsilon-phase  Hugoniot  of  RHA  are  listed  in  Table  1 and 
are  plotted  as  extensions  of  preceding  alpha-phase  Hugoniot  data  in 
Figures  6 and  7. 


IV.  DISCUSSION  AND  CONCLUSIONS 

The  Hugoniot  of  epsilon-phase  RHA,  shown  in  Figure  7,  is  displaced 
to  the  left  of  the  epsilon-phase  Hugoniot  of  iron  at  stresses  just  above 
the  phase  transition.  This  suggests  that  RHA  undergoes  a volume  offset 
at  the  phase  transition  which  slightly  exceeds  that  of  iron.  This  result 
is  consistent  with  results  reported  by  Takahashi  and  Bassett8  who  used 
X-ray  diffraction  at  static  high  pressures  to  show  that  the  volume  offset 
at  phase  transformations  in  iron-nickel  alloys  increases  with  nickel  con- 
tent (RHA  specimens  contained  from  3.0  to  3.5  percent  nickel).  At 
stresses  above  50  GPa,  the  Hugoniot  of  epsilon-phase  RHA  is  displaced 
above  the  Hugoniot  of  epsilon-phase  iron  indicating  stiffer  behavior. 

For  at  least  5.0  GPa  above  the  transition  stress  in  RHA  the  risetime 
of  the  P2  wave  is  approximately  0.3ns  based  on  the  5 to  95  percent 
criterion  adopted  by  Forbes3.  This  risetime  is  consistent  with  the 
corresponding  0.2  to  0.3us  risetime  of  the  P2  wave  in  iron  reported  by 
Novikov,  et.  al.9.  Barker  and  Hollenbach2,  and  Forbes,  and  it  suggests 
transformation  kinetics  similar  to  those  of  iron.  At  23.6  GPa  the  rise- 
time of  the  P2  wave  in  RHA  has  decreased  to  approximately  0.1ns,  again 
in  close  agreement  with  results  for  iron  reported  in  Reference  2. 

The  risetime  of  the  P2  wave  was  found  to  be  0.07ns  in  the  test  at 
34.5  GPa.  However,  this  risetime  may  not  be  accurate  because  the  PI 
and  P2  waves  appear  to  be  merging.  As  they  merge,  the  feature  identified 
as  the  top  of  the  PI  wave  may  be  displaced  to  higher  velocities  (See 
Figure  3).  In  this  test,  the  transition  stress  has  apparently  increased 
to  13.74  GPa.  Although  the  data  for  iron  by  Barker  and  Hollenbach  reveal 
an  increase  in  the  amplitude  of  the  PI  wave  with  increasing  driving  stress, 
it  can  be  shown  that  an  increase  to  13.74  GPa  is  unreasonably  large  for 
RHA  in  this  test.  Horie  and  Duvall30,  in  their  treatment  of  wave  propa- 
gation in  a phase  transforming  material,  assume  the  material  to  be 


8T.  Takahashi  and  W.  A.  Bassett , "The  Composition  of  the  Earth’s 
Interior ",  Scientific  American , Vol.  212,  No.  6,  June  1965,  pp.  100-8. 

9S.  A.  Novikov,  I.  I.  Divnov,  and  A.  G.  Ivanov,  " Investigation  of  the 
Structure  of  Compressive  Shock  Naves  in  Iron  and  Steel",  Sov.  Phys.  - 
JTEP,  Vol.  20,  No.  3,  March  1965,  pp.  545-6. 

10T.  Horie  and  G.  E.  Duvall,  "Shock  Waves  and  the  Kinetics  of  Solid- 
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PARTICLE  VELOCITY  (m/s) 


STRESS  (GPa) 


initially  overdriven  along  a metastable  extension  of  the  Hugoniot  for 
the  low-pressure  phase.  The  overdriven  material  then  relaxes  into  the 
high-pressure  phase.  In  the  Horie- Duvall  model,  stress  in  the  PI  wave, 
Op^,  changes  as  the  wave  travels  away  from  the  impacted  surface,  decreas- 
ing from  the  initial  overdriven  stress,  o,,  to  a stress,  a , , after  long 

a pia> 

travel  times.  Stress  in  the  PI  wave  is  given  by  the  expression. 


°pl  = °pl« 


+ (od  • aPi J exp  c-*/2*). 


where  t is  time  after  impact  and  t is  the  relaxation  time  which  is 
assumed  to  be  constant.  The  time,  t,  is  given  approximately  by  the 
expression, 

t = k CX1  - X0)/Upl.  (17) 

where  (Xj  - Xq)  is  the  specimen  thickness,  Up^  is  the  velocity  of  the  PI 

wave,  and  k is  a constant.  The  constant,  k,  appears  because  time,  t,  is 

less  than  the  transit  time  of  the  PI  wave;  the  reflected  elastic  wave  in 

RHA  reduces  the  stress  in  the  following  PI  wave  by  approximately  2.0  GPa 

and  no  further  attenuation  should  occur  after  o , < a , (See  Reference  2) . 

pi  pi®  J 

For  both  RHA  and  iron,  k ? 0.92.  Equation  16  represented  the  data  of 
Barker  and  Hollenbach  within  approximately  ±0.1  GPa.  However,  a t and 
°pl~  cou^  not  be  found  to  produce  this  level  of  agreement  with  the  data 

for  RHA,  and  the  problem  seemed  to  be  associated  with  the  location  of 
the  datum  point  of  Test  540  (ap2  = 34.5  GPa).  Since  the  risetime  of  the 

P2  waves  in  RHA  is  approximately  the  same  as  reported  for  iron,  it  was 
assumed  that  transformation  kinetics  might  be  similar  and  that  the  relax- 
ation time  for  iron  might  be  applicable  for  RHA.  The  results  for  RHA, 
using  t = 0. 17y s and  oploo  = 12.80  GPa,  are  shown  in  Figure  8.  The  data, 

except  for  the  datum  point  of  Test  540,  are  closely  represented  and  a 
stress  of  12.95  GPa  is  predicted  for  the  PI  wave  in  Test  540.  Data 
associated  with  this  adjusted  value  of  opl  are  shown  in  parentheses  in 

Table  1;  the  datum  point  for  the  epsilon  phase  undergoes  negligible  change 
as  a result  of  this  adjustment. 

The  epsilon-phase  Hugoniot  of  RHA  has  been  determined  at  stresses  up 
to  135  GPa.  There  is  close  agreement  with  the  Hugoniot  of  iron  over  the 
entire  stress  range,  with  only  minor  displacements.  The  volume  offset  at 
the  alpha-epsilon  phase  transition  in  RHA  is  slightly  greater,  and  the 
response  of  RHA  at  higher  stresses  is  slightly  stiffer.  On  the  basis 
of  risetime  in  the  P2  wave  and  relaxation  time  in  the  Horie-Duvall  model, 
the  transformation  kinetics  in  RHA  and  iron  should  be  similar.  Similar 
close  agreement  between  the  Hugoniots  of  RHA  and  iron  was  also  found  in 
the  preceding  study  of  the  alpha  phase  in  which  the  linear  (U-u)  relation- 
ships of  RHA  and  iron  were  essentially  parallel  and  displaced  predictably 
by  the  difference  in  bulk  velocities  at  room  temperature  and  pressure. 
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LIST  OF  SYMBOLS 
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t time,  ps  (subscripts  0,  1,  2,  3,  and  4 refer  to  times  in  Figure  2). 

u particle  velocity,  m/s  (subscripts  e,  pi,  and  p2  refer  to  the  elastic, 
PI,  and  P2  waves,  respectively. 

uQ  impact  velocity,  m/s. 

uf  free-surface  velocity,  m/s  (subscripts  1,  2,  and  4 refer  to  free- 
surfaces  velocities  in  Figure  2). 

U wave  velocity  with  respect  to  material  ahead  of  the  wave,  m/s 

(subscripts  e,  pi,  and  p2  refer  to  the  elastic,  PI,  and  F2  waves, 
respectively) . 

U*  wave  velocity  in  laboratory  coordinates,  m/s. 

3 

V specific  volume,  cm  /g. 

X distance,  mm  (subscripts  0,  1,  2,  3,  and  4 refer  to  distances  in 
Figure  3) . 

3 

p density,  kg/m  (subscripts  0,  e,  pi,  and  p2  return  to  unstressed 
material,  the  elastic,  the  PI  wave,  and  the  P2  wave,  respectively). 

o stress,  GPa  (subscripts  e,  pl,  and  p2  refer  to  the  elastic,  PI  and 
P2  waves,  respectively). 


Subscript  "i"  refers  in  a general  way  to  an  initial  condition. 
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